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We investigate the response of an electron system which exhibits ideal nesting features. Using 

the standard Matsubara formalism we derive analytic expressions for the imaginary and real parts 

of the bare particle-hole susceptibility. The imaginary part has sharp peaks whose maxima at the 

nesting momenta approximately scale with (lo/T). The peak lineshapes resemble neutron scattering 

data on chromium and some copper oxide superconductors. The real part of the bare susceptibility 

at the nesting vectors diverges logarithmically at low temperatures. Analytic formulas for the first 

O , vertex correction to the susceptibility are derived for a Hubbard interaction, and its momentum 

and temperature variations are calculated numerically. This term detracts substantially from the 

\-^ ' ordinary RPA terms for intermediate values of the Coulomb repulsion. Exact cancellation of a 

^—4- , certain class of diagrams at half filling is shown to result from particle-hole symmetry. We discuss 

the consequences of these results for spin fluctuation theories of high temperature superconductors 

and spin density wave instabilities. 
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I. MOTIVATION 



'"rt ! Nested Fermi surfaces with nearly parallel orbit segments are traditionally associated with electronic instabilities. A 

£J5 ■ peak in the susceptibility at a nesting momentum can create a charge density wave for strongpfilectron-phonon coupling, 

or alternately induce a spin density wave (SDW) when the Coulomb repulsion dominatesjj These instabilities often 

appear in systems with reduced-dimensionality for electron dynamics like the quasi one-dimensional organic metals 

and layered structure materials. a 

t-H , The present work on a square Fermi surface is motivated by the discoveries of nesting phenomena in high temperature 

superconductors. In these materials, the anomalous normal state quasiparticle damping that is linear in frequency and 

temperature can be explained by nesting, providing that electron-electron scattering is the primary damping source. 

The microscopic nested Fermi liquidu (NFL) theory derives this unconventional damping from electron scattering 

across nested regions of Fermi surface that approximately satisfy the nesting condition e(k + Q) + e(k) = 2/x, where 

Q is the nesting vector. This implies scaling of the spin susceptibility at the nesting vector in frequency divided 

[ — ' by temperature (ui/T) which has been observed by neutron scattering experiments. By contrast, the weak Fermi 

0^ ■ liquid damping that is the hallmark of ordinary metals owes its quadratic frequency and temperature variation to a 

^i ' susceptibility that is essentially independent of temperature. Within a self-consistent scheme, the NFL theory has 

provided a physical explanation of tlpp optical conductivity, electronic Raman spectrum, and marw other normal state 

features of cuprate superconductorsfl Fermi surface nesting has been discovered by photoemissiorJj in L^S^CaC^Os 

and some other cuprate compounds. 

Our present study is relevant to the spin fluctuation mechanism of d-wave superconductivity which was originally 
examined thirty years ago by Berk and Schrieffer.u The basic Fermi sphere yields extremely small transition tem- 
peratures T c for d-waves, and thus this concept became dormant until evidence for an unconventional mechanism in 
heavy fermion metals revived interest. The possibility of raising T c by.-a susceptibility enhancement treated in the 
random phase approximation (RPA) was suggested by Scalapino et alB and subsequently used by many groups in 
2D tight-binding model calculations foe cuprates. Resting enhances d-wave electron pairing in leading order exchange 
of antiferromagnetic spin fluctuations .0 Schriefferja however, has stressed the need to examine vertex corrections and 
other terms beyond the simple RPA. If the RPA series yields a large susceptibility enhancement, it is reasonable to 
expect that self energy and vertex contributions become corriespondingly important £9 

Within the nesting approximation, Virosztek and RuvaldaU have derived analytic expressions for the self energy 
and vertex corrections to the susceptibility at the nesting vector. These higher order contributions preserve the scaling 
of the susceptibility and they become comparable to the RPA terms for intermediate values of the Coulomb repulsion. 
Thus the RPA method by itself can be misleading, especially when the system is close to an SDW instability. 

In this paper we study these issues by using a 2D linearized dispersion model which gives a square Fermi surface at 
any filling. The linearization allows the imaginary and real parts of the noninteracting susceptibility to be obtained 
analytically and the perfect nesting allows us to examine the validity of the NFL theory. 



The linear dispersion model was first introduced by Mattial3 as a simplified-, version xif the 2D tight-hinding model. 
A similar model with square Fermi surface has also been studied by Luther £3 MattisEj and HlubinaEfl have claimed 
that the model is exactly solvable using bosonization and that the system is a 2D Luttinger liquid, with spin-charge 
separation, rather than a Fermi liquid. Indeed, bosonization is one of the powerful methods that can be used to solve 
the ID version of this model — which is a variant of the Tomonaga-Luttinger model.t^l In higher dimensions, however, 
bosonization (at least in the simple version used by Mattis and HAubina) is expected to give an incomplete picture 
of the system due to the presence of the particle-hole continuum.t3 HaldanelHl has pointed out that the distinctive 
feature of ID fermion systems is the nonexistence of low-energy particle-hole pairs with small momenta, whereas in 
higher dimensions this continuum always exists. This absence of excitation decay channels allows the response of a 
ID system to be completely described in terms of bosonic collective excitations. In the language of susceptibility, the 
imaginary part of the spin and charge susceptibilities in ID will consist only of delta function peaks corresponding 
to these bosonic modesJISI In higher dimensions, we get a qualitatively different picture. In addition to peaks from 
the collective excitations, there is also a regular contribution to the susceptibility, corresponding to excitations of 
particle-hole pairs which may also serve as a decay channel of the collective excitations if their dispersions lie within 
the particle-hole continuum. This feature appears clearly in our analytic results for the susceptibility, and may shed 
light to the current controversy concerning the existence of Luttinger liquids in dimensions higher than one.113 

We define the square model and our convention for units in Sec. II. In Sec. Ill we proceed to evaluate the imaginary 
and real parts of the bare susceptibility and display the key features emanating from this square model in the 
noninteracting limit. Effects of electron-electron interaction are studied in Sec. IV where we numerically compute the 
leading vertex correction to the susceptibility and compare it with the RPA term of the same order. In this section 
we also prove a particle-hole theorem that is relevant to a large class of diagrams near half filling. We discuss the 
results of our investigation in Sec. V. 

II. SQUARE FERMI SURFACE MODEL 

We consider a modified Hubbard model on a square lattice defined by the Hamiltonian 

h = h + Hu = 5>(k) - n) c + CT c k , CT + u J2 ni,-rma. (!) 

k,cr i 

with a noninteracting energy dispersion given by 

e(k)=v F (\k x \ + \k v \), (2) 

for momenta in the first Brillouin zone — Tr/a < k x ,k y < n/a. In the coordinate space the noninteracting Hamiltonian 
has the form 

oo 

^o = -E Yl *M ( c L c i +Xm , CT + 4,,c i+ym<a ), (3) 

i.er m=— oo 

with x m = (2m — 1, 0), y m = (0, 2m — 1), and hopping amplitudes 

t(m) = — £ -^ . (4) 

Note that the electrons are only allowed to hop from one sublattice to another, but not between sites in the same 
sublattice. 

In the rest of this paper we will use the following convention for units: vp = ir/a = kg = 1, where kg is the 
Boltzmann constant. The unit of energy is therefore vf(tt/o) = W/2, where W is the noninteracting bandwidth. If 
the bandwidth were 2 eV, temperature T = 0.03 would correspond to 360 K. The symbol /i will be reserved for the 
chemical potential. We also define a useful variable 7 = 1/(2T) = (3/2. 

The Fermi surface for this model with /1 = 0.9 is the square shown in Fig. |l|. The arrows specify the nesting vectors 
for this value of the chemical potential that will be used throughout this paper. 

Assuming the velocity to be independent of momentum avoids the logarithmic divergence in the density of states 
that characterize other square Fermi surface models (such as the 2D tight binding band at half filling) . This allows 
us to study the characteristic features due to nesting separately from those caused by the Van Hove singularities. For 
the present model, the density of states (per spin) is 



N(e) = e e(e)0(l - s) + (2 - e) 9(e - 1)6(2 - s). 



(5) 



We use the standard Matsubara finite temperature formalism to calculate the bare and the first vertex correction 
to the spin susceptibility of this model. The bare susceptibility is given by 



T 



Xo(q,m>) = -— ^2G (k,iv)G (k + q,iuj + iis), 



(G) 



k.iu 



where Go(k, iu>) = [iui — e(k)] 1 is the bare Green's function. Doing the frequency sum and performing analytic 
continuation to real frequency, we get the standard result for the bare particle-hole susceptibility 



,„ x 1 [^ /[g(k + q)]-/[ £ (k)] 

Xo(q,^ = - / dk — - — — — , 

4 J lu - e(k + q) + e(k) + id 



(7) 



where f(x) = [exp (3(x — ji) + 1] 1 is the Fermi function. At u> = and q — > this reduces to the temperature 
independent Pauli susceptibility x(0,0) = N(/j,). 



III. SUSCEPTIBILITY FOR NONINTERACTING ELECTRONS 



A. Imaginary Part of Susceptibility 

The susceptibility comprises real and imaginary parts xoio.,^) — Xo( f l' a; ) + ^Xoi^ 1 - )- The expression for the 
imaginary part is 



xo(q,w) 



dx 



dy{ 
i "• 



l 



e /3(s(x,y)-n) _|_ i e f3(e(x+ 



*,v+b)-») + 1 } 6 { W ~ £(x + a ' V + b) + £ ( X ^)}' ( 8 ) 



where a = max(|g x |, Iq^Q and fe = min(|g x |, \q v \) with < b < a < 1 for q G 1BZ. The integral can be carried out 
directly by dividing the region —1 < x,y < 1 into 16 regions separated by the lines x = { — 1, — a, 0,1 — a, 1} and 
y = { — 1, —b, 0, 1 — 6, 1}. The result is conveniently expressed in terms of the following functions: 



*(-)^{/m-/m}=-*(^ 



(9) 
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-S 3 (x,-u>), 



(11) 



K 1 (a,b,ui) = {a + b-\u)\)@{a + b-\u)\)@(\u>\-a + b) + 2b@(a-b-\u>\ 



(12) 



#2(0, b,ui) = @(a + b-u>) <3{u-a + b) + 0(a + b + u) Q(-u - a + b), 



(13) 



K 3 (a, b, w) = 0(a + b-u>) 0(uj + a - b) + 0(a + 6 + w) 0(-w + a - 6), 



(14) 



£>i 



(a, 6) = (^-) {s 2 (l - a) - 5 2 (1 + a) + S 2 (l - b) - S 2 (l + 6) + 5 2 (2) - S 2 (0) + S 2 (a + fe) - S 2 (2 -0-6)}, 



(15) 



Diia, b) = (^-) {s 2 (l + a - 6) - 5 2 (1 - a + 6) + S 2 (2 - a) - S 2 (2 - 6) + S 2 (6) - S 2 (o)}. 



(16) 



In the definition of S% (x) , we have used the dilogarithm function which is defined by 

Li 2 (a;) = dt ln(l - t)/t. (17) 

J X 

Using these functions, the imaginary part of the susceptibility can be written as 

Xd'(q^) = Xo,rcg(q, w ) +X0,sing(<L W )> (18) 

with a regular part 

Xo'rcg(q^) = {&(4 - a - b,w) + 5i(2 + a - &,w) + Si (2 - a + b,w) + S^a + 6,cj)} Ki(a,b,u) 

+ (s 3 (A - a- b,u)~ S 3 (2 + a-b,u) + S 3 {2 - a + b, u>) -S 3 (a + b, u)\ K 2 (a,b,u) (19) 

+ {#,(4 - a - 6, w) + £3(2 + a - b, w) - £3(2 - a + 6, w) - S 3 {a + &, w)} #3(0, 6, w) , 

and a singular contribution 

Xo,sing( c l' w ) = - D i( a ' 6 ) {^(w-o-6)-5(w + a+6)} + D 2 (a,b) {<J(w + a - 6) - tf(w - a + 6)}. (20) 

For small (lu/T) we can expand .5i(a;,a;) and S&(x, uj) 



Si(x,w) = (—J sech ^^^J + 0{{w/T)% (21) 

&(*,«) = (?£) [l + tanh(^^)+0(( W /T) 2 )]. (22) 

These expansions can be used to obtain the leading term in (w/T) of the regular part at the nesting vector Q = (/i, /i) 

W»"l " e < 2 « - M) (fi) { 1 + > -h'(^) h- ~*(^) + £ -K^) } + o(^) ■ W> 

For a half-filled system (/1 = 1), this reduces to 

«.«) = e( 2 -M)@+o@, (24) 

which is just the leading term of the exact expression that can be obtained using the nesting conditional 

x£ reg (Q,w) = (tt/2) JV(1 - w/2) tanh(a;/4T). (25) 

The momentum variation of %q (q, w) is sharply peaked at the nesting vectors as shown in Fig. g. Neu- 
tron scattering epqperiments have observed this type of peak structure in chromium and also in high temperature 
superconductors.fL2l The strong temperature variation of the peak lineshape is caused by nesting since the quasiparti- 
cle damping does not enter in the noninteracting susceptibility. A 3D plot shows the full momentum dependence of 
Xo,re g (q,w) in Fig. |. 

The imaginary part of the susceptibility at the nesting vector in Eqs. (p3|)-(p5|) approximately scales with frequency 
divided by temperature. The weak energy dependence of the density of states, which is nowhere singular in this model, 
introduces only a slow modulation which can be neglected for small oj, leaving only the temperature T tcij-set the 
effective energy scale. Scaling of the susceptibility was first obtained by means of the nesting approximation,Hu and it 
determines the anomalous quasiparticle damping that distinguishes the cuprates from ordinary metals. Fig. H shows 
the frequency variation of x'o reg(Q> w ) a ^ several temperatures. 

The singular part of the susceptibility comes from forward scattering processes which are strongly enhanced in this 
model by the synergetic combination of the linear dispersion and the flatness of the Fermi surface. These terms are 
key ingredients for collective charge and spin modes and they are relevant to the issue of charge and spin separation. 
In Fig. [5] we show the frequency dependence of Xo (q> w ) f° r small q. The delta function peak carries a spectral weight 
which is two orders of magnitude larger than the integrated weight of the continuum. Nonlinear energy dispersion 
or finite curvature of tho-Fermi surface, however, will remove this singularity and replace it with a very sharp, but 
nonetheless finite, peakEJ The 2D nearest-neighbor tight-binding model, for example, gives a finite susceptibility for 
all q except for the half- filled case when there is a Van Hove singularity in the density of states. 



B. Real Part of Susceptibility 

The real part of the susceptibility determines the SDW instability and is a vital component in the exchange processes 
that may create a d-wave superconducting state for a system with highly anisotropic Fermi surface. Taking the real 
part of Eq. (M) we obtain the expression for the real part of susceptibility 

, (l\ f 1 f 1 tanh[ 7 (e{x, y) - n)] - tanh[ 7 (e(x + a,y + b)- /i)] 

Xo(q>w) = (-) V / dx dy — ■ r . (26) 

v8/ y_ x y_j u-e{x + a,y + b) + e{x,y) 

Both integrals are principal value integrals as indicated by the symbol V '. The integral can be done by dividing the 
region of integration into 16 regions as was previously done for the imaginary part. The result can be expressed in 
terms of functions R\ and R 2 which are defined in the Appendix. We hrst define several auxiliary functions: 

L\{a, b) = ln[cosh7(2 — a — fi)] — ln[cosh7(2 — /j,)] — ln[cosh7(2 — a — b — /j,)] + ln[cosh7(2 — b — fi)] 

+ ln[cosh7(l — fi)] — ln[cosh7(l — b — fi)] + ln[cosh7(l + a — b — /j,)] — ln[cosh7(l + a — jj)] 

+ ln[cosh7(l — fi)] — ln[cosh7(l + b — //)] + ln[cosh7(l — a + b — /i)] — ln[cosh7(l — a — fi)] 

+ ln[cosh7(a — /i)] — ln[cosh-f(a + b — /.i)] +\n[cosh-f(b — /j,)] — ln[cosh.7(— [/,)], 



(27) 



L 2 {a,b,uj) = (—) [i?i(l; 1-a + b; 2 - a + b - to) + i?i(l + a - b\ 1; 2 + a-fr-w)} 
+ (— ) |i?i(a; b; a + b - w) + Ri{2 - b: 2 - a; 4-a-b-u/)}, 

L 3 (a,b,u>) = ( a+ ~^ ) {^i(2-a; 2-a-b: 4 - a - b - u>) + R x (2\ 2-6; A-a-b + uo)\ 

[i?i(l-a + 6; 1-a; 2 - a + 6 - w) + Ri{l + b; 1; 2-a + 6 + w)| 

[i?i(l + a; l + a-6; 2 + o-6 + w)+i?i(l; 1-6; 2 + o-6-w)| 
( a - j [i?i(6; 0; a + 6-w) + i?i(a + 6; a; a + 6 + w)}, 



(28) 



32 

a + 6- 


a; 


32 

a + b- 


UI 


32 

a + 6- 


U! 



(29) 



4-0-6-w) + i? 2 (2; 2-6 
2-a + 6-w) + ^2(1 + 6; 1 
2 + a-6 + cj) - R 2 (l; 1-6 



Li(a,b,uj) = i?2(2 — a; 2 — a — 6 
+ i? 2 (l - a + 6; 1-a 

- R 2 {l + a; 1 + a-b 

— R 2 {b] 0; a + b — uj) — R 2 (a + b; a; a + b + uj) 
With these functions, the real part of the bare particle-hole susceptibility is 

'D 1 (a,b)\( 1 1 i (D 2 {a,b) 



4- a -b + uo) 
2-a + b + uj) 
2 + a-b-u)) 



(30) 



KCU,) = (r/4)L l( a,6) + (^^){^^--^^-} + (^^)(^^ L^\ 

v// v; V 7T /lw + a + 6 w-a-6J V tt / I w - a + 6 w + a-6J 

+ \ L 2 (a,b,Lu) + L 2 (a,b,-Lu) + L 3 (a,b,uj) + L 3 (a,6,— w) + (T/4) (L 4 (a,6,w) + i 4 (6,a,tj))L 



(31) 



The momentum variation of the real part of the susceptibility exhibits broad peaks at the nesting vectors, as seen 
in Fig. pi for a frequency uj = 1 meV and bandwidth W = 2 eV. At room temperature these nesting peaks exceed the 
Pauli susceptibility by roughly a factor of two. 

For a = b = fi = 1 (nesting vector of a half-filled system) the expression for the real part of the susceptibility 
simplifies to 

X' Q (Q, w) M =i = -2Tln[cosh(l/2T)] + (1^)^(1,0, 2 - u) + (^±^) J R 1 (1, 0, 2 + to). (32) 

At the static limit and T — ► this asymptotically goes to 

x[ ) (Q,0) M =i = -2Tln[cosh(^)]+ln(|^), (33) 

where je = 1.781... is the Euler constant. The susceptibility at these limits can also be derived more easily using the 
density of states and the nesting condition, which is exactly satisfied at half-filled.u 

The logarithmic divergence-that appears in the susceptibility is vital to the SDW and the competing d-wave super- 
conducting state instabilitiesEil that need to be carefully examined when we turn on the electron-electron interaction. 
In the next section we proceed to the study of vertex corrections to the susceptibility. 



X rpa (q> *9o) = ^ TTZJZ ;„_\ - ( 34 ) 



IV. COULOMB INTERACTIONS 

In this section we will use the standard diagrammatic analysis to discuss the effects of the interaction Hamiltonian 
Hjj in Eq. (|J) . The Hubbard onsite Coulomb repulsion U only couples electrons with opposite spin on the same lattice 
site. This, in particular, means that all exchange diagrams vanish since there is no interaction between electrons with 
equal spin. 

A. Random Phase Approximation 

The RPA is one of the most useful approximation methods in many-body physics and has been used effectively to 
describe plasma oscillations and density wave instabilities. In this approximation, the longitudinal spin susceptibility 
for the interacting Hamiltonian is given by the simple formula 

2xo(q i ?: go) 

1 -Uxo(<l,iqo)' 

The diagrams corresponding to the first three terms in this series are shown in Fig. OTa)-(c) where the bare suscepti- 
bility "bubble" is (a). Note that for each diagram there is another diagram in which all the spins are flipped. In the 
paramagnetic state, which will always be assumed in the following calculations, this spin degeneracy gives rise to the 
factor of two in Eq. (|3J) . 

Despite its popularity, the RPA should be used with caution since occasionally its application can lead to serious 
mistakes whose remedy requires a study of vertex corrections. For example, one may incorrectly infer that a single 
component electron gas should be superconducting when the screened Coulomb interaction is treated by this method, 
but Sham et al. have shown that .the true situation of a net repulsion between electrons is confirmed when vertex 
corrections are taken into account E3 |— ■ 

For our discussion, it is useful to define the spin-dependent particle-hole susceptibilitiesE3 

X ffff '(q,^o) = ^ J dr e^ T (T T [n CT (q,r)n CT ,(-q,0)]), (35) 

where 

n ff (q,0) = ^ C j c+qa c kiCT . (36) 

k 

From this susceptibility we can get the charge and spin susceptibilities 

x c (q%) = x TT (q%) + x u (q%) + 2x n (q,%), (37) 

x s (q,%) = x TT (q,%) + x u (q,%) - 2x U (q,«9o)- (38) 

At zeroth order in U, we have x^(q^9o) = X^(q>*9o) = Xo{<i,iqo) and X^H^^Qo) = 0. It is useful to define separate 
RPA series for these susceptibilities 



,TT („ in \ - Xo(q^go) _ V U 

l-(C/Xo(q%)) 2 



Xr T pa(q> i( lo) = : _ (TJ ', . ,^ 2 = Xrpa(q> «?(>), (39) 



y n (q ian) = ^(xo(q^o)) 2 iT ( . x f401 

Within RPA, the charge and longitudinal spin susceptibilities are therefore 

Xrpa(q^%)) = 2(x r T p r a (q,i9o) + X r T pa(q>i<?0)), (41) 

Xr pa (q, m) = 2(x r T p r a (q, m) - x r T pa(q> m))- (42) 

The RPA would give an exact susceptibility if in the formula we replace the bare susceptibility Xo(qi* ( 7o) with 
the exact irreducible susceptibility which requires all possible self energy and vertex corrections. Clearly this is very 
difficult. For the tight-binding model, the spin and charge susceptibilities have been calculated to second order in 
U.t3 The calculation, however, was done on a finite lattice and the susceptibilities were calculated only at q = (0, 0) 
and q = (it /a, ir/a). In the next subsection we will discuss the momentum and temperature dependence of the leading 
order vertex correction to the bare particle-hole susceptibility for the square model. 



B. Vertex Correction to X (l>*9o) 
The first vertex correction to x (Piipo) 1S shown in Fig. H(d) and is given by the expression 

Xcorr(p,ipo) =-"AfX,X, r (KP,iko,m) G (k + p,ik + ip Q ) G (k,ik ), (43) 



N 

k ikn 



where 



U 2 T 
r(k,p,ifc ,ip ) = -jj-^2^2xo{q,m) G (k + q,ik +iq ) G (k + p + q, ik + ip + iq Q ). (44) 

q iqo 
Note that this is the only vertex correction to X (li *9o) hi the second order. To calculate this correction we define 

U 2 

t . 

k 

with 



Xcorr(p,«Po) = -— ^F(k,p,ip ), (45) 



i /013 r< 2 3 r 114 _i_ r< 24 

^-.-w-^s V-Vte-w - (46) 

G q/3 _ T 2 y^ 1 y^ xo(q, »?o) , 47 % 

where a € {1,2}; /? e {3,4}; £1 = £ k+q ; £ 2 = (£k+ P + q - ipo); 6 = £k; & = (£k+ P - ipo)- Note that g and p are 
boson frequencies, k$ is a fermion frequency, and Xo(q, *<7o) = Xo(l) w ~* *9o)- The two frequency summations in G a ^ 
can be performed using the analytic expressions for Xo(q, *9o) that we have obtained previously. First note that for 
a simple pole l/(iqo — z) we have 



t 2 y — - — y - 



, iko - 6<9 • iqo + iko - Za iqo - z 

lk iq ^ 

[coth( 7 z) - tanh( 7 ^)] p h (7(^ Z *» Z^tM }, 

Now observe that in the expressions for xo(q, *9o) each term has exactly one simple pole of the form l/(iqo — Zi) where 
Zi varies from term to term. (Li(a, b) and £3(0, b, iqo) clearly are not in this form, but one can check easily that their 
sum can be written in this form.) Except for the poles multiplying Di(a,b) and D2(a,b), all the other poles reside 
in an integral. We assume that the integration and the frequency summation operations commute so that we can 
perform the frequency summation first. With all this, G a/3 can be obtained from Xo(Q, iqo) by a simple substitution 
rule 

[coth(7Zi) - tanh(7£ a )] 



i<7n — z i 

(49) 
r tanh(7(g a - z t )) - tanh(7Jg) \ 

l £,a - Zi - £/3 V 

Performing this substitution, G a @ is now expressed as a sum of one dimensional integrals. These integrals can 
be calculated analytically using the piecewise polynomial approximation to the functions tanh(a;), ln[cosh(x)], and 
coth(x) as presented in the Appendix. Note that the number of terms that need to be calculated in this method 
increases very rapidly compared to the number of terms in the bare susceptibility. The analytic expressions, however, 
allow us to perform the principal value integrations exactly within this approximation. With this analytic expression 
for G Q ^, the computation of Xcorr(p, ipo) is reduced to a numerical computation of a four dimensional integral. 

The integral over q was carried out numerically to give the function F(k, p,0) defined in Eq. (ffq). We have 
performed these computation using SP2 parallel processors at the Maui High Performance Computing Center and the 



University of Virginia. Each of the plots of -F(k, p, 0) shown in Fig. || typically requires about 4 days of CPU time 
on 10 processors. The vertex correction Xcorr(P)O) is then obtained by integrating F(k, p,0) over k. 

In Fig. I) we show the result for Xcorr(p, 0) at T = 0.03 at several points in the first Brillouin zone. The influence of 
this vertex correction to the static susceptibility is demonstrated in Fig. y by the dashed curve, in comparison to the 
bare susceptibility (dotted curve) and the second order RPA term (dot-dashed curve). Combining these corrections 
yields the solid curve for the total susceptibility which exceeds the bare curve by roughly a factor of two when U = 0.5. 

The temperature variation of the Xcorr(p,0) at the nesting vector Q = (n, p.) is displayed in Fig. fLO. The bare 
susceptibility (solid curve) and the second order RPA (dot-dashed curve) diverge at low temperatures as ln(T) and 
In 3 (T) respectively. The standard RPA spin susceptibility Eq. ( p4| ) yields a Neel temperature for an SDW instability 
at Tjv = 60 K for W = 2 eV and U = 0.5 eV. The vertex correction Xcorr(p, 0) is expected to diverge also as ln 3 (T) 
and we have fitted the computed points (circles) using a third order polynomial in ln(T). 

C. Vertex Corrections to x (q>*9o) 

At second order in U, there are two vertex correction diagrams to x (l>*?o) which are shown in Fig. [111. The 
expressions for these diagrams are 



\A 



(9) = C/2 EEE G °( fc i) G "( fc i + «) g o(p) Go(k 2 -ki+p) G (fc 2 ) G (fc 2 + q), (50) 



fei fc 2 



Xa(q) = ^ 2 EEE G °( fc i) G °( fc i + «) °o(p) G o( fc 2 + h - p) G (fc 2 ) G (k 2 - q) (51) 

fci fc 2 P 

where we have used q = (q, iqo) and 

q iq a q 

These diagrams are usually expressed in terms of the bare particle-hole susceptibility Eq. (o) and the particle-particle 
susceptibility 

<Ao (q, *w) = jf ^2 G °( k ' ^) G o(q - k, fw - iv), (53) 

We remark that at half filling these two susceptibilities are related by 

0o(Q-q,«*>)=xo(q,«w), (54) 

where Q is the nesting vector at half filling: Q — (1,1). For the following discussion we will use the forms given in 
Eq. (|50|)-(|5l|). The fact that the Hubbard interaction vertex does not depend on momentum directly implies 

XA(q,%) =Xcorr(q«Qo), (55) 

therefore our previous results for Xcorr(q, *9o) ar e still useful for this section. The other diagram XB(q, iqo) differs 
from XA(q, *<?o) only in that the particle lines in one of its two loops are replaced by hole propagators. In general these 
two diagrams must be calculated separately. We show, however, that at half filling, where we have exact particle-hole 
symmetry, these two diagrams exactly cancel each other. To show this, we write 

XA(q) + XB(<Z) = U 2 J2J2 G °( fcl ) G °( fcl + «) G °(P) A ( fcl <^' «)' ( 56 ) 

fel V 

A(ki,p,q) = J2k 2 G o(h-k 1 +p)Go(k 2 )G (k 2 + q) 
+ E fc2 Go(fc 2 + fci - p) G (k 2 ) G (k 2 - q). 

The symmetry between the two terms can be seen more clearly if in the second sum we write k 2 = —k-j 



A(fci,p,g) = Ek 2 Go(+(k2-h+p)) G (+fc 2 ) G (+(fc 2 + g)) 
+ E k3 Go(-(h - fci + P )) Go(-h) G (-(*3 + g)). 

The momenta in the second sum are just the negative of the momenta that appear in the first sum. Now consider a 
propagator G(—k) = (—iko — £(— k)) _1 = (— iko — ^(k)) _1 . If we perform the following transformation 

iko^ipo] k^p + Q; f(k)->-£(p); (59) 

we get G(—k) — * (~ipo + £(p)) _1 = —G(p). This means that if we do the same transformation in the second sum 
in Eq. (pq) by writing ik^o = ik^o and k3 = k4 + Q, each propagator becomes the negative of the corresponding 
propagator in the first sum. Since there are three propagators in each sum, this directly implies that 

A(ki,p,q)=0. (60) 

We thus conclude that xa(q, iqo) + Xb(*1, iqo) — for a half-filled system. This cancellation is valid for any external 
momentum and frequency. For non half-filled systems, we expect the sum to vanish in powers of (/i/W) where p, 
is the chemical potential measured from the center of the band, and W is the bandwidth. The absolute scale is 
determined by xa(q, *9o) which, as we have shown, is exactly equal to Xcorr(q, «9o) which has been analyzed in the 
previous section. 

D. Particle- hole Symmetry 

The cancellation of diagrams presented in the previous section clearly can be generalized to higher order diagrams. 
This leads to a simple theorem: for a Hubbard interaction at half filling, a diagram which has a loop consisting of 
n > 1 particle lines is equal to (—1)" times a similar diagram where the direction of all particle lines in the loop is 
reversed. There are three things to note: (a) The theorem does not hold for the simple Hartree loop (n = 1) but 
applies to any other particle loop, (b) This cancellation depends only on the particle-hole symmetry which is exactly 
satisfied at half filling. It is therefore valid for any symmetric dispersion in any dimension, (c) In any diagram, particle 
lines which are not continuously connected to the external lines must reside in a loop. Therefore the theorem has 
relevance for a large class of diagrams. 

As a simple application of this theorem, consider the self-energy diagrams in Fig. n2. The second order diagram in 



Fig. 12(a) is the first nontrivial self-energy diagram in a Hubbard model. The expression for this diagram is 

X(k) = U 2 J2Go(k-q) X o(q)- (61) 
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For a nested Fermi surface, Virosztek and Ruvalds have shown that this diagram yields a quasiparticle damping which 
is linear in frequency and temperatureE The RPA method, which is usually used to selectively include higher order 
diagrams, replace the bare susceptibility Xo(qO with Xrpa(<?)- The self-energy in this approximation is given by 

oo 

£ rpa (fc) = ££gS>(£0, (62) 

ra=l 

Eg:'(fc) = t/ 2 » ]T G (k - q)xl n - 1 (q)- (63) 

i 

Note that there is no odd order term in this series. The fourth order RPA self-energy, T, rp ' a (k), is shown in Fig. |l2](b). 
For a Hubbard interaction, which does not depend on momentum, this diagram is exactly equal to another self-energy 
diagram due to the particle- hole ladder with four rungs, T. h (k), which is shown in Fig. |l2|(c). In general, we have 

zg:Hk) = x^\k). (64) 

At half filling, we can use the previous theorem to relate the self-energy diagrams due to the particle-hole ladder 
with the corresponding diagrams from the particle-particle ladder. The exact relation is 

4h ) ( fc ) = (- 1 ) m 4p ) ( fc )- (65) 



Combining these results, we conclude that at half filling 

Xgf(k) = Z < £ ) (k) = Z% n \k), (66) 

£jr +1) (fc) + 4 2 p " +1) (fc) = 0. (67) 

Using these relations, together with £r P 2 = 0, we can write the sum of the contributions to the self-energy from 
the RPA, particle-hole ladder, and particle-particle ladder at half filling to be 

£ sum (fc)=[/ 2 ^G (fc-q) X sum(<z), ( 68 ) 

q 
where 



3xo (g) 

1 - (Uxo(q)) 2 



Xsum(q) = 1 ,™: H ' - 2xo(<z). (69) 



We have substracted 2xo(?) since there is only one diagram for the three series at the second order. 

V. DISCUSSION 

The present results for a square Fermi surface in general support the nested Fermi liquid theory of the anomalous 
damping that distinguishes high temperature superconductors from ordinary metals. Analytic forms for the imaginary 
part of the susceptibility reveal a sharp peak structure in momentum space which will assure the dominance of electron 
collisions across nested regions. Scaling of the susceptibility as a function of lo/T is modified slightly by the energy 
dependence of the density of states, but otherwise fits the requirements to produce a linear frequency and temperature 
variation of the quasiparticle damping. 

The real part of the susceptibility for the square model demonstrates a broader peak shape in momentum space and 
its maxima at the nesting vectors display a slow divergence in ln(T) that is important for the spin density wave and 
the competing d-wave superconducting phase instabilities. The analytic results presented in this paper may provide 
a basis for estimating the transition temperatures for Hubbard models. 

The cumbersome task of evaluating the multidimensional integrals for the first vertex correction is slightly simplified 
by the analytic formulas of the bare susceptibility and some approximations that we develop for the required hyperbolic 
functions. The numerical evaluation of the resulting expressions, however, still requires considerable computing power 
and time. The calculated vertex correction counteracts the susceptibility enhancement that comes from the simple 
RPA and therefore should not be neglected in computations of d-wave superconducting transition temperatures. 

Particle-hole symmetry in a half-filled energy band leads to a theorem for the cancellation of a general class of 
diagrams for a Hubbard interaction. This theorem eliminates the need to evaluate diagrams containing loops with an 
odd number of propagators at half filling. 

In contrast to the standard BCS theory in simple metals, where Migdal's theorem allows the neglect of certain vertex 
corrections for the electron-phonon coupling, spin fluctuation theories normally cannot rely on a disparity between the 
energy scales for the interaction medium and the bandwidth. Hence the reliability of many calculations for the d-wave 
superconducting transition temperatures are difficult to estimate because vertex corrections for tight-binding models 
are particularly difficult to compute. Our proof for cancellation of many types of higher order diagrams applies also 
to the tight-binding models at half filling and it is reasonable to expect that these diagrams will be less important 
even in the realistic cases such as the copper oxide superconductors. Nevertheless it would be interesting to derive 
these corrections in terms of the chemical potential placement away from half filling. 

A quantitative constraint on our weak-coupling analysis is the strength of the Coulomb repulsion U which should 
remain much less than the bandwidth. For £//W <C 1, the combined RPA and first vertex correction yield a net 
enhancement of the susceptibility. Higher order vertex correction terms clearly should be included for larger U. Self 
energy corrections have been examined using the nesting approximation in the nested Fermi liquid theoryJ3 The self 
energy corrections for the square model will be explored in future work. 

Finally, the controversial issue of charge and spin separation in two dimensional systems with nesting may benefit 
from the complete analytic results for the bare susceptibility which includes the regular part in addition to the singular 
terms emanating from the flat sides of the Fermi surface. Both terms should be included in future computations of 
the electron self energy, with particular emphasis on particle-hole decay channels and their impact on the spectral 
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function. Similarly these susceptibility formulas may feature new charge and spin collective modes which may be 
relevant to a variety of interesting metals. 
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APPENDIX 

The functions R\ and i? 2 arc defined by (the dependence on 7 and fj, are implicit) 

R t {a- b- c) = M^ia - l x), 7 (c/2 - /,)) - M^(b - „), 7 (c/2 - ft)), 1 g {1, 2}, (70) 

with 



Jo z-w 



(71) 



,. l x lncosh(z) , 
M 2 (x, w)=V I — dz. 



(72) 



These integrals cannot be expressed in terms of known functions. In this paper we have chosen to analyze these 
integrals by approximating the functions tanh(x) and ln[cosh(x)] 



tanh(x) — > x9(ci — |x|) + \ [ 1 + m(\x\ — c 2 ) 3 ] 6(|x| — ci)6(c 2 — \x\) + Q(\x\ — C2) > sign(x), 



(73) 



ln[cosh(x)] - — 9(ci - \x\) + {\x\ + -(\x\ - c 2 ) 4 - c 3 } 9(\x\ - Ci)0(c 2 - \x\) + {\x\ - c 3 }9(|x| - c 2 ), (74) 

with c\ = (4/25); m = (625/11907); c 2 = (67/25); C3 = (1691/2500). Using these approximations, the integrals can 
be analytically calculated. For Mi(x,w) we obtain 



Mi{x, w) = Q(x, w) Q(x) - Q(-x, -w) Q(-x), 



(75) 



Q(x, w) = Qi(x, w) @(ci — x) + Q 2 (x, w) 9(x - ci)0(c 2 — x) + Qz{x, w) Q(x — c 2 ), 



(76) 



3 l 

)(x,w) = y^-(io - c 2 ) 3_ "(x - c 2 ) n + (w-c 2 ) 3 ln\x -w\, 



n=l 



(77) 



Qi(x,w) = x + w In 



w 



i(x,w) = Qi(ci,w) +m (Q (x,w) - Qo{c\,w)) + In 



C\ — w 



Q 3 (x, w) = Q 2 (c 2 , w) + In 



c 2 - w 



Similarly, for M 2 (x,w) we get 



M 2 (x, w) = B(x, w) Q(x) + B(-x, -w) ©(-a;), 



(78) 
(79) 
(80) 

(81) 
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B{x, w) = B\(x, w) 0(ci — x) + B2(x,w) Q(x — Ci)@(c2 — x) + B^(x, w) Q(x — c 2 ), 



4 1 
B (x,w) = V -(«, - c 2 ) 4 -"(x - ca) n + (w - c 2 ) 4 In |.t - ■ 



(82) 
(83) 



a; — w 






B 2 (x,w) = Bi(a,w) + (to/4) (S (a;,w) - S (ci,w)) + (a; -ci) + (w - c 3 ) In 



Ci — 10 



(84) 
(85) 



S 3 (a;, to) = S 2 (c 2 , to) + (a; - c 2 ) + (w - c 3 ) ln 



c 2 - 10 



(86) 



Finally, in the calculation of the vertex correction, we also need to approximate the function coth(x). We do this 
using the following polynomial 



coth(x) -> (l/x)6(c 4 - \x\) + sign(x)(l + c 5 /x 2 )0(\x\ - c 4 ), 
with C4 = 1/2 and C5 = 1/4. 



(87) 
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FIG. 1. Fermi surface for the square model using a chemical potential fi 
this Fermi surface. 



0.9. The arrows denote the nesting vectors for 
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FIG. 2. Xo,rc g ( c l) w ) shows sharp peaks at the nesting vectors Qi = (fi, fi) and Q2 = (2 — /j,, 2 — /u) at low temperatures. The 
model parameters are chosen to be W = 2 eV, w = 1 meV, and /1 = 0.9 eV. The three curves are for 2T/W = 0.01, 0.02, and 
0.03. These nesting peaks contrast with isotropic electronic structure results which yield a susceptibility that has a smooth 
momentum dependence. 

FIG. 3. The momentum dependence of Xo,reg(k,u>) for [i = 0.9, u> = 0.01 and T = 0.03. The peaks at the nesting vectors are 
strongly temperature dependent. 

FIG. 4. Frequency variation of Xo,reg(Qi w) is shown at 60 K (solid curve), 120 K (dashed curve), and 240 K (dotted curve) 
using W = 2 eV and fi = 0.9 eV. 

FIG. 5. Frequency dependence of Xo(Qj w ) f° r Q = (.QiQ) consists of a regular particle-hole continuum and a singular peak. 
The peak has the form A(q)6(w — 2q) with A(q) increases linearly with q as shown in the inset. 

FIG. 6. Real part of the susceptibility x'(k, w) in the static limit (u> = 0) as a function of momentum shows the nesting 
peaks. We use T — 0.03 and \x = 0.9. Note that the vertical scale begins from 0.9, which is the density of states at the chemical 
potential N(fi). 

FIG. 7. The first three terms of the RPA series (a-c) and the second order vertex correction diagram (d) for x TT ■ Note that 
diagram (b) belongs to x 1 ^- 

FIG. 8. Clockwise from top left: the function F(k,p, 0) plotted as a function of k for p = (0.1, 0), (0.5, 0), (0.9, 0.9), and 
(0.9, 0). 

FIG. 9. The static susceptibility as a function of momentum showing the bare susceptibility x(q>0) (dotted curve), second 
order RPA contribution U 2 x 3 (q, 0) (dashed curve) , vertex correction Xcorr (q, 0) , and the total of these three terms (solid curve) . 
U = 0.5 (W/2). 

FIG. 10. The static susceptibility at a nesting vector x(Q>0) diverges as ln(T) at low temperatures as shown by the solid 
curve. The second order RPA contribution f/ 2 X 3 (Q> 0), which diverges as ln 3 (T), is shown by the dot-dashed curve. The circles 
represent the results of our numerical computation for Xcorr(Q,0). This vertex correction is expected to diverge as ln 3 (T) and 
we have drawn a fit (dotted curve) for the data using a polynomial 3.766 + 3.483a.- + 0.9530a; 2 +0.08844a a where x = ln(2T/W). 
We have used /j, = 0.9 (W/2) and U = 0.5 (W/2). 

FIG. 11. The first two vertex correction diagrams for X (Qi *8o)- The diagrams only differ in the orientation of lines in the 
second loop. These two diagrams exactly cancel each other at half filling. 

FIG. 12. Diagrams (a) and (b) are the first two self energy diagrams in the RPA; diagram (c) is the third term in the self 
energy series due to the particle-hole ladder; diagram (d) is the third term in the particle-particle ladder. For a Hubbard onsitc 
interaction diagrams (b) and (c) are equal and at half filling diagram (d) is also equal to diagram (c). 
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